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Pixel correlation searches for O VI in the Lyman a forest 
and the volume filling factor of metals in the Intergalactic 
Medium at z ~ 2 — 3.5 
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ABSTRACT 

Artificial absorption spectra are used to test a variety of instrumental and physical 
effects on the the pixel correlation technique for the detection of weak O VI absorption. 
At H I optical depths ^ 0.3 — 1, the apparent O vi detections are spurious coincidences 
due to Hi absorption at other redshifts. In this range, the apparent O vi optical depth 
is independent of Hi optical depth. At larger Hi optical depths, apparent Ovi optical 
depth and H I optical depth are correlated. Detailed modelling is required in order 
to interpret the significance of this relation. High-resolution spectra of four QSOs 
together with a large suit of synthetic spectra are used to show that the detection of 
Ovi in individual spectra is only statistically significant for overdensities >5. These 
overdensities are larger than would be naively inferred from the onset of the correlation 
and a tight optical depth-density relation. The lower limit for the volume filling factor 
of regions which are enriched by O vi is 4% at 95% confidence. This is no larger than 
the observed volume filling factor of the winds from Lyman break galaxies. Previous 
claims that the observed O vi absorption extends to underdense regions and requires 
a universal metal enrichment with large volume filling factor, as may be expected from 
population III star formation at very high redshift, appear not to be warranted. 

Key words: intergalactic medium - quasars - galaxies:formation - metal enrichment. 



1 INTRODUCTION 

The Lycv forest in QSO absorption spectra is now widely 
believed to be caused by fluctuating Gunn-Peterson absorp- 
tion due to an undulating warm photoionised Intergalactic 
Medium (IGM). With this interpretation for the origin of 
the Lya forest, it has become possible to relate the absorp- 
tion optical depth to the density of absorption systems in 
a meaningful manner. According to this scheme low column 
density absorption systems arise from low density regions of 
the Universe. A detection of associated metal absorption in 
weak Lya absorption systems would suggest that the metal 
enrichment of the IGM is widespread. 

Models put forward to explain this metal en- 
richment fall into two categories: they either propose 
that early and widespread Population III star forma- 
tion in pre-galactic structures at z ^ 10 is responsible 
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or suggest a later episode of metal enrich ment by winds 
from starbursting galaxies at z < 5 (e.g. lAguirre et alJ 
1200 ll iTheuns et alTl200ll lAdelberger et alJ l2003ft .' Trans- 
porting metals from the galaxies or protogalaxies where 
they are produced into the (low density) IGM is not, 
however, a trivial matter. A variety of mechanisms 
have be en proposed. Som e have argued for supern ova- 
driven dCouchman fc Reesl Il98rj. iDekel fc Silkl Il98rj) or 
'miniquasar '-driven iHaiman. Madau fc Loeblll999l) winds. 
Others have argued fo r ej ection via g alactic mergers 
jGnedin fc Ostrikerlll997l and|Gnedin|| l998|) or ph otoevapo- 
ration during reionisation iBarkana fc Loeblll999f) . With an 
accurate determination of the fraction of the universe (or 
'volume filling factor') enriched by metals it may be pos- 
sible to test the veracity of such models. A large volume 
filling factor would favour early enrichment by stars form- 
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ing in shallow potential wells, out of which metals are more 
easily transported into low density regions. 

Observationally, the most important tracers of met- 
als in the photoionised IGM are C IV and O VI . Associated 
C IV absorption has be en detected down to HI column densi- 
ties of order 10 14 cm" 2 llTvtler et alll99llCowie et alll995l 
lEllison et ai]l2000l iPettini et alJl200St corresponding to an 
optical depth of 1 — 10. To investigate even lower density 
systems a pixel correlation technique has been employed 
where, instead of fitting an associated C IV absorption line, 
a pixel-by-pixel search for excess absor ption at the corre- 
spond ing C iv wavelength is performed (ICowie fc Songailal 
1998). This method was hoped to be more robust than 
the u se of line 'stacking' llTvtler et al.lll995l | Elliso n et alJ 
2000), but there is some ambiguity in the interpretation of 
the results. At the relevant redshifts, the O VI species is the 
most sensitive tracer of metals in low density regions of the 
IGM, provided th at the spectrum of the UV b ackground 
is sufficiently hard (iRauch. Haehnelt fc Steinmet3ll997l and 



Table 1. Model parameters used for all simulations 



iHellsten et alJll998f) . However, the relevant Ovi absorption 
lines (Xovm = 1032A, A vit = 1038A) occur in the same 
wavelength range as the Lyman a forest and higher order 
lines of the Lyman series, which makes the detection of weak 
Ovi absorption c h alleng ing. 

ISchave et alJ i2000f) claim to have detected Ovi at 
Hi optical depths significantly smaller than unity in some 
QSO absorption spectra and used the optical depth-density 
relation of numerical hydrodynamic simulations to argue 
that this indicates a detection of metals in underdense re- 
gions of the IGM. In this work, we apply the pixel- by- 
pixel search to synthetic spectra and thereby investigate in- 
strumental and physical effects on this technique (see also 
lAguirre. Schave fc Theunsl l2002f) . We then assess the sta- 
tistical significance of the Ovi detection and the implied 
density threshold above which the presence of metals is con- 
fidently detected. 

The spectra (both simulated and observed) are de- 
scribed in Section|21 In Section|3]we outline the pixel correla- 
tion method used to detect the presence of Ovi absorption. 
The results of the search for Ovi are shown and discussed 
in Section 01 In Section |S] we explore the consequences for 
the metal enrichment of the low density IGM. 



2 SIMULATED AND OBSERVED SPECTRA 
2.1 Calculating synthetic spectra 

2.1.1 Basic Assumptions 

We use the method develo ped by Bi tei. Boerner fc Chul 
Il992l: iBi fc DavidsenI Il997l) to analytically calculate syn- 
thetic spectra in the fluctuating Gunn-Peterson approxima- 
tion for a warm photoionised IGM with a density distribu- 
tion that is expected in a ACDM structure formation model. 
The basic assumptions are: 

(i) the baryonic matter traces the dark matter on scales 
larger than a filtering scale, which is related to the Jeans 
scale; 

(ii) the distribution of mildly non-linear densities has a 
lognormal probability distribution function (PDF); 

(iii) the absorbing gas obeys a simple temperature- 
density relation and is not shock-heated. 



Parameter Chosen Value 



h 

log[T eB \ 
V 



0.3 

0.65 
0.7 
0.9 
4.6 
0.641 
1.333 



Table Q shows the model parameters. A brief summary 
of how we calculate the synthetic spectra follows. 

2.1.2 The Spatial Distribution of Baryons 

According to standard structure formation models the mat- 
ter distribution was initially a Gaussian random field with 
small density fluctuations that evolve under gravity into the 
pattern of sheets and filaments characteristic of CDM mod- 
els. Many features of such a distribution can be reproduced 
by a rank-ordered mapping of a linear Gaussian random field 
onto a lognormal PDF of the density. This idea is at the 
heart of the method us ed to calculate synthetic absorption 
spectra iBi et alJ il992). This method is an efficient way of 
producing large numbers of realistic synthetic spectra. Note 
that we are interested in producing synthetic spectra with 
a realistic density PDF and realistic instrumental proper- 
ties. Reproducing the detailed clustering properties of the 
density field is less important. 

We begin by creating a linear Gaussian random 
field which is fully determined by its power spectrum, 
PpM(k), for which we have u sed the form given by 
Efsta thiou. Bond fc White! \l2 



i(*) 



(1) 



(1 + [ak + (Wc) 3 / 2 + (cfc) 2 ]") 2 /" ' 

where a = 6.4/r, b = 3.0/r, c = 1.7/r, F = Q. m h (Q m 
& h have their normal definitions) and v — 1.13. The power 
spectrum was normalised by the rms fluctuation amplitude 
on a 8h _1 Mpc scale, as, as given in Table 1. In order to take 
into account pressure effects that suppress fluctuations on 
scale s smaller than a certain filtering scale jGnedin fc Huil 
1998) the power spectrum of the baryon density in the linear 
regime is assumed to have the form, 



Ps(fc) 



2 7 fcT cff 



3/im p f2(l + z) 



(2) 



(l + a; 2 fc 2 ) 2 ' 

where Xb is a comoving scale re lated to the Jeans length 
(Peebles 1980; Bi & Davidson 1997), 7 is the ratio of specific 
heats, T e g is the effective mean temperature and is the 
molecular weight of the IGM. An appropriate value for T e s 
has been deter mined from compariso n with hydrodynamic 
simulations bv lBi fc DavidsenI <ll997f) . This is the baryonic 
matter power spectrum at z = and linear theory is used 
to evolve it backward to high redshift. 

Equations Q and © provide the 3D power spectrum. 
A ID power spectrum needed for the description of line-of- 
sight (LOS) fluctuations is obtained by an integration of the 
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Table 2. Observed Data Sample 



QSO 


z cm 


z rangc 


Dovi 


D„, 


Q1122-165 


2.40 


2.02-2.34 


0.157 


0.165 


Q1442+293 


2.67 


2.51-2.63 


0.194 


0.193 


Q1107+485 


3.00 


2.71-2.95 


0.254 


0.253 


Q1422+231 


3.62 


3.22-3.53 


0.418 


0.386 



3D power spectrum jKaiser fc Peacocklll99ll) . The matter 
density and the peculiar velocity field are related and can 
be de scribed as coupled Gaussian random fields fsee lBi et alJ 
Il992l for more details). Random realisations of the density 
and peculiar velocity in real space (co-moving coordinates) 
are obtained using a Fourier transform routine. The non- 
linear evolution is modelled by a local rank-ordered mapping 
to a lognormal density probability distribution, 



n/n — exp 



(3) 



where 8 — (n — n)/n and n/n is the overdensity fac- 
tor. At the relevant redshift the baryonic density field is 
still in the mildly non-linear regime and t he use of a log- 
normal PDF is a reaso nable assumption (iBi et al.l Il992t 
iNusser fc Haehneltll2000l) . 



2.1.3 From the density to the absorption spectrum 

The Hi optical depth depends not only on the density 
but also on the temperature of the gas, due to the tem- 
perature dependence of the recombination rate. For den- 
sities n/n < 10, numerical simulations show that most of 
the gas is not shocked and that a simple power law re- 
lation between density and temperature is established by 
the balance of photoioni sation heating and adiabatic cool- 
ing jHui fc Gnedinlll997l) . 



T = T {n/n) c 



(4) 



where 0.3 < a < 0.6 and T is a constant. 

For gas at temperatures of a few times W 4 K in pho- 
toionisation equilibrium the neutral hydrogen density is pro- 
portional to n 2 T _0 ' 7 /rjj j, where Tjjj is the photoionisa- 
tion rate. As a result, the optical depth to Lya absorption 
at redshift z can be written as 



component of the Hubble parameter and oo is the cross- 
section for resonant Lya scattering. 

The effect of peculiar velocity and Doppler broadening 
are taken into account to obtain the optical depth in redshift 
space. We have thereby approximated the Voigt profile by a 
Gaussian distribution, which is a reasonable approximation 
for absorption lines with equivalent width < 0.7A. 

The optical depth for the O VI doublet is calculated as- 
suming a fixed ratio of O VI number density to H I number 
density (novi/'iHi)- The optical depth of the higher order 
Lyman lines are also calculated. 

2.1.4 Realistic Synthetic Spectra 

To facilitate a more accurate comparison with the observed 
spectra we model the optical depth distribution for the same 
wavelength range as the observed spectra. There are two re- 
gions of the spectra of particular interest; the 'Ovi region', 
where Ovi absorption is searched for, and the 'Hi region', 
where the corresponding Lya absorption occurs. The spec- 
tra are obtained by co-adding the optical depth of Lya, 
Ovi and higher order Lyman series lines. We then use the 
estimated statistical error on the flux for the observed spec- 
trum to add random noise and perform a Gaussian smooth- 
ing to mimic the effect of instrumental broadening. 

The optical depth is scaled such that the mean flux 
decrement of the observed QSO is reproduced. As we will 
discuss later the results of the pixel-by-pixel search depend 
sensitively on the mean flux level. We have thus rescaled 
the optical depth for the Hi and the Ovi regions indepen- 
dently to reproduce the observed mean flux level in both 
regions {D m and D \i respectively). D m is used to set TLya 
in the Hi region as well as fovi (with fixed novi/fim) and 
the optical depth of the higher order Lyman lines in the 
O vi region. -Dovi is then used to set fL ya in the O VI region. 
Note that this means that we allow a difference in the value 
of fi^ar/rHi in the two regions of the spectrum. It seems 
plausible to assume that the errors we obtain in this way 
for a range of Monte Carlo realisations of the LOS density 
distribution are similar to those we would have found if we 
had only picked random realisations that can reproduce the 
mean flux decrement in the H I and the O VI regions simulta- 
neously. The latter approach is computationally prohibitive. 

In this way we have obtained ensembles of simulated 
spectra with the same wavelength range, mean flux decre- 
ment and noise properties as the observed spectra, but dif- 
fering in their specific random realisation of LOS matter 
distribution and noise. 



where 



I{z) = 



cdz 



TL ya (z) = a nml{z), 



1 1 \ \ 1.7 
n{z ) N 



1 1+2 HoE(z') V n 



■ exp 



(•») 



(6) 



z\ and Z2 are the lower and upper limits of the redshift 
region of interest, b is the Doppler parameter, nm is the 
neutral hydrogen density, v pcc is the peculiar velocity and 
a value a = 0.4 is assumed. E(z) is the redshift dependent 



2.2 Observations 

We have used four observed QSOs which are part of the 
sample used bv lSchave et al.l (|200CT) for a detailed compar- 
ison with synthetic QSOs spectra: Q1122-165, Q1442+293, 
Q1107+485 and Q1422+231. The data for Q1122-165 was 
taken during Commissioning I and Science Verification ob- 
servations for the UVES instrument on the VLT (Kueyen) 
and released by ESQ for public use. The reduction m ethod 
can be found in iKim. Cristiani fc D'Odoricol (1200 ll) . The 
other spectra were t aken with the HIRES instrument 
iVost si I S. et al.lll994]}_ on Keck I and r educe d using pro- 
cedures described in Barlow fc Sargentl lll997l) . These two 
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Figure 1. An example of the results from the pixel-by- pixel 
search method for the absorption spectrum of Q1107+485. 



instruments provide spectra with resolutions -R max ~ 80000 
and 110000 respectively and have typical S/N of 50. Table 
[5] gives a list of the redshift range and flux decrement for 
the Ovi and Hi regions. Regions of strong absorption from 
other known metal lines were removed. 



3 PIXEL-BY-PIXEL SEARCH 

We use a technique for the de tection of Ovi in the Lym an 
a forest initially developed bv lCowi e fc Smigailj|j[l99 8|) for 
the detection of weak Civ absorption. ISchave et alJ (2000) 
adap ted this metho d to search for the O VI doublet (see also 
lAguirre et alJ [20021 . We briefly summarise the procedure 
here. The redshift range explored in each of t he four QSOs 
is show n in Tabled This is the same as used bv lSchave et alJ 
( 2000). Due to the increase in noise we did not find it worth- 
while to extend the search to lower wavelength. 

The continuum-fitted absorption spectrum is re-binned 
to pixels of equal size in redshift space for the Lyman se- 
ries and both members of the Ovi doublet. A width of 
(cAA/AL y ri)kms -1 j s taken (where Lyn is the highest or- 
der Lyman line used). In order to take into account possible 
velocity shifts between O VI and H I absorption, we have also 
tried larger pixels widths. This did not improve the sensi- 
tivity for O VI detection. 

A search for correlated absorption at the wavelengths 
of the Lyman series and the O VI doublet is then performed. 
Each pixel of the Lya absorption region is considered in 
turn. If the flux level is within <7 n oi/2 of the continuum level 
(where <r no i is the estimated statistical error on the flux for 
the observed spectrum) or in a region where the associated 
Ovi absorption is not covered the pixel is discarded. If the 
pixel is within <T no i/2 of saturation we attempt to use the 
higher order Lyman lines to estimate TL ya . We thereby use 
the higher order line with the lowest equivalent Lya optical 
depth TL yQ = min(rLyn/LyQ//LynALyn) that is not within 
Cnoi/2 of saturation or the continuum (where fhyn is the 
oscillator strength of the Lyn line). 

Once we have established a value for TLy a , we consider 
the pixels corresponding to absorption by the O VI doublet 
at the same redshift. The limiting factor in the use of 
this method is coincident absorption from Hi, either in 



the form of higher order Lyman lines at similar redshift or 
by Lya lines at lower redshift. Obviously we want to min- 
imise this contamination. The doublet nature of the O VI is 
utilised in this regard. The smaller equivalent absorption, 



tqvi = mm rovta 



/OviaAoviaTOvi 
/ovibAovib 



(7) 



will be the least contaminated. This gives us an up- 
per limit to the O VI optical depth. We will thus call 
the Ovi optical depth obtained in this way the apparent 
O vi optical depth. 

As pointed out by lAguirre et alJ i2002T) the contam- 
ination by higher order Lyman lines can be estimated 
from the corresponding Lya absorption and corrected. Such 
contamination is subtracted from the O VI optical depth 
where possible. If this is not possible because of saturated 
Lya absorption we also discard the pixel. 

We then bin the ensemble of optical depth pairs in 
Lya optical depth and plot the relation between the median 
optical depth of O VI and Lya . We use the median instead 
of the mean as the distribution of O VI optical depth in each 
bin is skewed toward high optical depths due to the contam- 
ination by H I absorption. 

Fig. Qshows an example of the relation of the apparent 
O VI and Lya optical depth (rovi,a PP — TL yQ relation) ob- 
tained in this way. At small Lya optical depth, the median 
Ovi optical depth is constant. As we will demon s trate in 
detail in the next section (see also lAguirre et all (|2002)), 
the apparent O VI is here almost exclusively due to contam- 
ination of coincident Hi absorption and this region of the 
diagram cannot be used to infer a detection of O VI . At 
larger Lya optical depth, the apparent Ovi optical depth 
rises and is mainly due to the presence of Ovi . The level 
of metal enrichment is, however, not straightforward and 
requires detailed modelling. 



4 RESULTS FROM SIMULATED SPECTRA 
4.1 Principle parameters 

To obtain an indication of how the pixel-by-pixel search de- 
pends upon instrumental effects and physical parameters 
and what errors should be expected due to variation of 
the LOS density distributions, we investigated a large suit 
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Figure 3. The rovi.app — Tjya relation for an ensemble of forty synthetic spectra. Left: Dovi increases from bottom to top. Right: Djji 
increases from bottom to top 



of synthetic spectra. We have varied the mean flux decre- 
ment in the Ovi and Hi regions (Do vl and -Dm), the ratio 
of O vi to Hi density (novi/nm) and introduced a varying 
cut-off in the density below which we set the O VI density 
to zero ((n/n) C ut )■ Fig. |21 sketches how the relation of ap- 
parent Ovi optical depth and Lya optical depth depends 
on the three key parameters Dovi, {n/n) cu t and novi/fim- 
A detailed discussion of the effect of these and the other pa- 
rameters using synthetic spectra of Q1122 will follow below. 
The parameters were varied around fiducial values denoted 
by an 'F' in figurcs lHlol The error bars in these figures are the 
la spread of the optical depths for ensembles of 40 spectra. 

4.2 Varying the mean optical depth 

Fig. |2l shows how changing the mean flux decrement af- 
fects the rovi.app — Thya relation. In the left panels we con- 
sider changes in the Ovi region, Do VI . With increasing 
flux decrement the level of spurious coincidences due to 
Lya absorption from gas at a lower redshift rises and mim- 
ics the detection of Ovi . Since this Lya absorption is un- 
correlated with absorption in the Hi region, this leads to 
a floor of constant apparent O VI optical depth which rises 
with increasing Do VI . In the right panels of Fig. |2]we vary 
the mean flux decrement in the Hi region Dm- This time 
the shape of the relation between H I optical depth and ap- 



parent O VI optical depth does not change but the error bars 
increase substantially with decreasing mean flux decrement. 
This is because, for small mean flux decrement, large optical 
depth Lya pixels are poorly sampled. 



4.3 Varying the O VI distribution 

In Fig. 0] we vary the O VI distribution for our synthetic 
spectra. The left panels shows how varying the ratio of 
O VI to Hi density affects the rovi,a PP — i~Lya relation. With 
increasing ratio a correlation between apparent O VI optical 
depth and Hi optical depth develops for large Hi optical 
depth but the errors are large. As discussed in the introduc- 
tion the main aim of the search for O VI at small H I optical 
depth is to establish the volume filling factor of metals in the 
Universe. Because of the possible spurious detections due to 
random coincidences this is not straightforward. 

There is another important effect which can lead to spu- 
rious detection of O VI . Low Lya optical depth is not neces- 
sarily associated with low density gas but can also be due to 
the wing of a strong Lya absorption feature which is caused 
by a high density region. In order to test the density range 
for which the presence of O VI is required to reproduce the 
observed rovi,a PP — T ~Lya relation, we introduce a cut-off in 
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Figure 4. As Figure 131 Le/t: novi/nni increases from top to bottom. Right: (n/n) cu t increases from top to bottom. 



the density, (n/n) cut below which we set the O VI density to 
zero, i.e. 



novi „ ., n n 
= 0, if - < - 

uhi n Vn/cut 



novi , ., n (n 

= const, it — ^ — 

uhi n \ n / cu t 



(8) 



(9) 



In the right panel of Fig.|I]we vary (n/n) cu t. The onset 
of the correlation moves to larger Hi optical depth. Unfor- 
tunately, this dependence is weak. 



4.4 Varying S/N and continuum level 

The left panels of Fig. |S]show the effect of changing the sim- 
ulated noise. As expected, the errors increase significantly 
with increasing noise level. 

The effect of an error in the continuum fit is shown in 
the right panels of Fig. where we have raised and lowered 
the continuum as a whole by 1% in each synthetic spectrum. 
When we raise the continuum level the apparent O VI optical 
depth increases because the optical depth of spuriously coin- 
cident H I absorption increases. Note that if the continuum is 
placed too low the correlation between Hi and Ovi optical 
depth misleadingly appears to extend to smaller H I optical 
depth. 
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Figure 5. As Figure 13 Left; noise increases from top to bottom. Right: the continuum level increases from top to bottom. 



5 SEARCH FOR O vi IN THE LOW DENSITY 
IGM 

5.1 Comparison of real and simulated spectra 

Here we perform a detailed comparison between four ob- 
served QSO spectra and sets of synthetic spectra. We use 
the same S/N, Dovi, Dm and wavelength range as in the 
observed spectra with one exception. In the case of Q1122 
we decrease Dovi by 35% in order to reproduce the level 
of apparent rovi at low thi- The need for this is most likely 
due to a deficit of saturated regions in our synthetic spec- 
tra. The contribution of saturated regions to the mean flux 
decrement is high at low redshift and the failure to ade- 
quately reproduce the high incidence of these systems leads 
to a bias in the mean flux decrement JViel et alJl200l . 

The triangles in the left panel of Fig. [§] show the 
Tovi.app — TLy Q relation for the four observed spectra. For 
the synthetic spectra we assume as before that the 
O VI distribution takes the form described in equations l|8)l 
and © . We have varied the O VI density and the overden- 
sity threshold for addition of O VI over the range novi/^Hi = 
— 0.3 and (n/n) cu t = — 100. We have produced samples 
of 40 synthetic spectra for each set of values. We have then 
calculated the rovi.app — TLya relation for each sample and 
assessed the agreement with the observed relation by cal- 
culating x 2 - Note, however, that there may be additional 



sources of error such inhomogeneities in the metal distribu- 
tion which our Monte Carlo technique does not take into 
account. The crosses and squares in the left panels of Fig.[S] 
show the rovi,ap P — TLya relation with no O VI and the best 
fitting values for novi/^m ((n/n) cu t = 0), respectively. The 
best fitting value has a reduced Xi close to or somewhat 
smaller than one. There is thus good agreement between 
the Tovi.app — T Lya relation for our best fitting simulations 
and the real data. In the right panels of Fig. |S] we use the 
Tovi.app — T ~Lya relation of synthetic and observed spectra 
for Q1107 to provide an example of the difficulty in con- 
straining (n/n)cut • 

The top panel of Fig.|7|shows the reduced x? as a func- 
tion of novi/nm for the four observed QSOs. In Q1442 there 
is a marginal detection of Ovi, in Q1107 O VI is detected 
with a poorly constrained Ovi density (wovi/whi ~ 0.08) 
while in Q1122 Ovi is clearly detected with novi/rom ~ 
0.06. In Q1422 no Ovi is detected. Indeed, we find that the 
spectrum of Q1422 is inconsistent with O VI absorption at 
the same level as detected in Q1122 and Q1107 with a con- 
fidence of greater than 99%. 

In the bottom panel of Fig. [7|we show how the reduced 
X 2 varies as a function of (n/n) cu tWe have assumed that 
novi/nm and (n/n) cu t are independent in the calculation of 
the confidence level. (n/n) cu t of greater than 4 for Q1122, 
7 for Q1107 and 4 for Q1442 are ruled out with 95% confi- 
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Figure 6. Right: the rovi,app — T^ya relation for observed and simulated synthetic spectra for four QSOs; this is shown for synthetic 
spectra with both the best best fitting constant ratio of no V i/"m a- n d with no O VI . Left: the effect of varying (n/n) cu t on synthetic 
spectra of Q1107 compared to the Tovi.app — Tjya relation of the observed spectrum. 
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Figure 8. The volume filling factor (top panel) and the mass 
fraction (bottom panel) as a function of gas overdensity for a 
lognormal density PDF at the redshift of the four observed QSOs. 



umn density Hi absorption systems in the spectrum of 
Q1122 again consistent with our results. 



dence. There is thus no significant detection of O VI at over- 
densities ^ 5. 



5.2 Other searches for Ovi absorption in the 
observed sample 

ISchave et alJ ll200(J) have claimed that the correlation of 
apparent Ovi optical depth and Hi extends to thi ~ 0.1 
in the case of Q1122, Q1107 and Q1442. They have in- 
terpreted this as a detection of O vi in underdense re- 
gions with n/n ~ 0.3 — 0.5. We cannot confir m this 
with our more detailed analysis. iDave et alJ l)l998i) found 
that Q1422 is consistent with no Ovi in a search for 
weak O VI absorption lines consistent with our result. 
Note that they detected Oviin a Ovi-C IV pixel corre- 
lation search which pro bed large Lyq optical depth. In 
the spectrum of Q1442, ISimcoe. Sargent, fc Rauchl l|2002l) 
found two strong Ovi systems associated with high col- 
umn density H I absorption systems but no weak ab- 
sorption systems consist e nt wi th our marginal detection. 
ICarswell. Schave fc Kiml l|2002l) found a large number of 
Ovi absorption systems associated with low and high col- 



5.3 Implications for the metal enrichment of the 
IGM 

The lack of a significant detection of O vi at low densities 
cannot be used as evidence against the presence of these 
metals. As discussed extensively in the previous section, 
Ovi at these densities may be masked by Hi absorption. 
Also, at high redshift the spectrum of the UV background 
may not be hard enough to ionise oxygen up to O vi. The 
decreasing contamination by Hi absorption (due to an ex- 
pected lower TLyof) is likely to be the main reason why the 
lowest redshift QSO in the sample is the only one for which 
we clearly detect Ovi at moderate overdensities. 

Fig. |Sl shows how the volume filling factor depends 
on overdensity for the lognormal density PDF in the red- 
shift ranges of the observed QSO absorption. For this 
PDF the 95% confidence limits for the density threshold, 
(n/n) cut , translate into lower limits of 4%, 1.5% and 4% for 
the volume filling factor of Ovi for Q1122 (z = 2.0 - 2.3), 
Q1107 (z = 2.7 -3.0) and Q1442 (z = 2.5 - 2.6), respec- 
tively. There appears to be no evidence for a large volume 
filling factor of metals from the observed oxygen absorp- 
tion. A lower limit of 1.5 — 4% for the volume filling factor is 
no larger than that inferred for winds from Lyman break 
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galaxies dAdelbereer et al.ll20od) . A picture where metal 
enrichment is due to winds from rather large galaxies at 
z ~ 2 — 5 is co nsistent wi t h the obse rved Ovi absorption 
in QSO spectra. lHaehnerd (ll99Sl) and iPettini et all <2003ft 
find that the same is t r ue for the observed Civ absorption, 
although ISchave et alJ i2003tl may find new evidence of rel- 
evance. 

Fig.|H|also shows the relation between mass fraction and 
overdensity for the lognormal density distribution. Despite 
the rather low inferred volume filling factor for which metals 
have been detected, about 20-40 % of baryons must already 
be enriched with metals to explain the observed absorption. 



6 CONCLUSIONS 

We have used a large sample of synthetic spectra which 
mimic the observational properties of four observed QSO 
spectra to interpret the results of a search for O VI in the 
low density 1GM. Our results can be summarised as follows. 

(f) At low Lya optical depth the corresponding appar- 
ent O VI optical depth obtai ned by using the pixel c orre- 
lation technique de veloped bv lCowie fc Songaila! (Il998h and 
ISchave et alJ a200Cf) depends mainly on the mean flux decre- 
ment in the O VI region of the spectrum. 

(2) In a sample of four QSO in the redshift range 2-3.5 
we detect O VI significantly in two QSOs and marginally in 
the third. The significance of the detection increases with 
decreasing redshift. 

(3) The position of the bend in the relation of apparent 
O VI and Lya optical depth depends only weakly on the low- 
est density at which OVI is present in the IGM. 

(4) We obtain upper limits of 4, 7 and 4 for the min- 
imum density for which Ovi has been detected with 95% 
confidence. For the lognormal model density distribution 
this translates into lower limits of 4%, 1.5% and 4% for 
the volume filling factor of metals. We thus do not confirm 
previous claims of a detection of O vi in underdense regions 
and a corresponding large volume filling factor of metals. 

(5) The O VI absorption in QSO absorption spectra as de- 
tected by the pixel correlation technique provides no ev- 
idence for (or against) a widespread metal enrichment at 
very high redshift (z ~ 10 — 20). The lower limit for the vol- 
ume filling factor of metals is equally consistent with metal 
enrichment by winds from Lyman break galaxies. 
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